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摘 要 : 为 提高 非 均匀 噪声 下 波 达 方向 〈direction ofarrival，DOA) 角 估 计算 法 的 估计 精度 和 分 辨 率 ， 基 于 低 秩 矩阵 恢 
复 理论 ， 提 出 了 一 种 二 阶 统计 量 域 下 的 加 权 L1 稀 政 重 构 DOA 估计 算法 。 该 算法 基于 低 秩 矩阵 恢复 方法 ， 引 入 弹性 正 
则 化 因子 将 接收 信号 协 方差 矩阵 重 构 问 题 转换 为 可 获得 高 效 求 解 的 半 定 规划 〈semidefinite programming，SDP) 问题 以 
重 构 无 噪声 协 方差 矩阵 ; 而 后 在 二 阶 统计 量 域 下 利用 稀疏 重 构 加 权 工 1 范 数 实现 DOA 参数 估计 。 数 值 仿真 表明 ， 与 传 
统 MUSIC、L1-SVD 及 加 权 Ll 算法 相 比 ， 所 提 工 法 能 显著 抑制 非 均匀 噪声 影响 ， 具 有 较 好 的 DOA 参数 估计 性 能 ， 且 
在 低 信 骂 比 条 件 下 ， 所 提 算 法 具有 较 高 的 角度 分 辨 力 和 估计 精度 
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Low-rank matrix recovery based DOA sparse reconstruction method 
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Abstract: To increase the accuracy and resolution of DOA estimation algorithm in the presence of non-uniform noise, based on 


the low-rank matrix recovery theory, this paper developed a weighted L1 sparse reconstruction DOA estimation algorithm (LR- 


WLOSRSS) under the second-order statistical domain. Following the low-rank matrix recovery approach, the algorithm firstly 


introduced the elastic regularization factor to transform the covariance matrix reconstruction problem of the received signal into 


a semidefinite programming (SDP) one that can be solved very efficiently to reconstruct the noise-free covariance matrix. After 
that, the DOA can be complemented in the second-order statistical domain by employing the sparse reconstruction weighted L1 
norm. Numerical simulations show that, compared to the traditional MUSIC, L1-SVD (Ll-norm-singular value decomposition) 
and weighted L1 methods, the proposed algorithm can suppress the influence of the non-uniform noise significantly, provide 
better DOA parameter estimation performance, as well as improve angle estimation accuracy and resolution in the case of the 
low SNR (signal-noise ratio) . 
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子 空间 类 算法 的 性 能 将 会 急剧 下 降 甚 至 失效 。 为 了 进一步 提高 

DOA 的 估计 精度 和 分 辨 力 ， 基 于 压缩 感知 理论 ，Liang G 等 人 
在 雷达 ,移动 通信 无线 传 感 器 网 络 和 射电 天 文学 等 领域 ， 利用 目标 信号 的 空域 稀疏 性 ， 提 出 一 种 稀 玻 重 构 1-SVD 算法 四。 
波 达 方向 角 估 计 (direction of arrival, DOA) 是 一 类 重要 的 问题 ， 在 信 源 数量 已 知 条 件 下 ， 该 算法 将 得 到 DOA 高 精度 估计 。 然 
通常 被 称 为 定向 (direction finding, DF) 估计 或 DOA 估计 口 。 而 ， 在 非 均匀 高 斯 噪声 或 未 知 信 源 先 验 信息 条 件 下 ， 该 算法 将 
基于 信号 和 噪声 子 空间 的 传统 DOA 估计 算法 显著 改善 了 DOA ”无 法 分 辨 两 个 邻近 和 角度， 因而 空间 分 辨 力 较 差 。 
的 估计 精度 和 分 辨 力 ， 如 多 重信 号 分 类 (multiple signal 近年 来 ， 随 着 信号 处 理 技术 的 快速 发 展 ， 众 多 研究 人 员 相 
classification, MUSIC ) 算法 让。 以 MUSIC 方法 为 代表 的 子 空间 类 提出 非 均匀 高 斯 噪声 下 的 DOA 估计 算法 。 文 献 [5] 提 出 一 种 
类 算法 通常 假设 噪声 是 零 均 值 、 方 差 为 1 的 复 高 斯 白 噪声 。 然 。”” 非 均匀 高 斯 噪声 下 最 大 似 然 (maximum likelihood, ML ) 估计 算 
而 当 附 加 在 阵列 传感器 上 的 噪声 为 非 均匀 ee 对 接 法， 该 算法 通过 逐步 迭代 求解 信号 和 噪声 的 对 数 似 然 函数 以 实 
收 信号 协 方差 进行 特征 空间 分 解 会 引起 信号 子 空间 泄露 ， 诸 如 。” 现 DOA 参数 估计 。 然 而 ， 众 所 周知 的 ML 算法 具有 的 较 强 初 
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录用 稿 


始 值 依赖 性 及 较 大 运算 复杂 度 限 制 该 算法 的 广泛 应 用 。 文献 [6] 
提出 一 种 非 均 匀 高 斯 噪声 下 基于 和 矩阵 补 全 的 DOA 参数 估计 算 
法 (matrix completion based music, MC-MUSIC)。 该 算法 基于 和 矩 
阵 补 全 思想 重 构 出 无 噪声 信号 协 方差 , 而 后 采用 传统 类 MUSIC 
算法 实现 DOA 参数 估计 ， 不 仅 降低 了 非 均匀 噪声 的 影响 ， 也 
进一步 避免 DOA 估计 算法 的 迭代 求解 。 然 而 ， 该 算法 没有 考 
虑 协 方差 矩阵 元 素 之 间 相 关 性 ， 有 可 能 导致 算法 求解 数值 不 稳 
定 ， 从 而 使 得 算法 稳定 性 较 差 。 

针对 上 述 问题 ， 本 文 基于 低 秩 矩阵 恢复 5 理论， 提出 一 种 
二 阶 统 计量 域 下 的 加 权 7 稀 琉 重 构 DOA 估计 算法 (Low-rank 
matrix recovery based Weighted J/ Sparse Reconstruction DOA 


Wr 


estimation algorithm under the Second-order Statistical domain, 
LR-WLOSRSS )。 非 均匀 高 斯 噪声 下 , 本 文 首 先 基 于 低 秩 矩阵 恢 
复 理 论 ， 引 入 弹性 正则 化 因子 将 接收 信号 协 方差 矩阵 重 构 问 题 
转换 为 可 获得 高 效 求解 的 SDP 问题 以 重 构 无 噪声 协 方差 矩阵 ， 
从 而 避免 非 均匀 高 斯 噪声 的 影响 ， 且 提高 了 低 秩 矩阵 恢复 算法 
求解 的 数值 稳定 性 ;而 后 在 二 阶 统 计量 域 下 利用 加 权 7 范 数 
(Weighted 1-norm ，WLI) 外 实现 DOA 参数 估计 。 数 值 仿真 
表明 , 与 传统 MUSIC、1-SVD 子 空间 类 算法 及 稀 疏 重 构 加 权 
算法 相 比 ， 所 提 算 法 能 显著 抑制 非 均匀 噪声 影响 ， 具 有 较 好 的 
DOA 参数 估计 性 能 ， 且 在 低 信 品 比 条 件 下 ,所 提 算 法 具有 更 高 
的 角度 分 辨 力 及 更 优 的 角度 估计 性 能 。 


1 ”接收 信和 号 模型 
假设 有 O 个 远 场 窄带 信号 {s (Dj2, 入 射 至 具有 jy 个 阵 元 
的 均匀 线性 阵列 ， 则 阵列 接收 信号 模型 可 表示 为 


z(]= 立 afe ja (+a 人 0) 
xG) 为 Mxl 信 号 接收 矢量 ， 


其 中 : $s(1?) 和 


a(0,)=[1 e* ev 分 别 为 第 9 个 信号 源 的 幅度 和 阵 


列 导 向 矢量 ，w = 2xsin(0,) /4 表示 阵 元 之 间 的 第 gq 个 信号 到 达 
此 阵 元 时 的 相 移 ，{6.,0,…,0,,…,6。} 为 @ 个 窄带 信号 的 发 射 角 
度 ，4 为 阵 元 间距 ，4 为 载波 信号 波长 ， 通 常 4 < 22， 

nD) =[n(D),w(D)…,nw (DJ 为 互 不 相关 的 0 均值 非 均 匀 高 斯 噪 
声 ， 即 n(1)~ CN(0,WW)，W 为 噪声 mb 功率 协 方差 矩阵 中， 信 


号 s,(1) 互 不 相关 。 
为 便于 推导 ， 式 (1) 接 收 信号 模型 可 进一步 重 写 为 
x(D)=As(D) tna(n) 多 


其 中 ， 阵 列 导向 流 型 矩阵 4=|e(Q),a(&)…a(bo)]sCe ， 


且 M > 8@， 即 假设 阵 元 数量 远大 于 信号 数量 ， 信 号 波形 矢量 


s(1)= [s (站 ,5 (站 so (HF ce Cox 。 


多 快 拍 数 条 件 下 ， 式 (2) 可 进一步 表示 为 
YL ~ 4VMx2S92xL + NY (3) 
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其 中 ;X=[x(0),x(2)…x(D]、 $=[s(),s(2).…s(D] 


N =[a(D,ma(2)…:a( 万 ] 分 别 为 工 个 快 拍 下 的 接收 信号 和 矩阵 、 


信号 幅度 矩阵 、 非 均匀 高 斯 噪声 矩阵 ， 世 表示 快 拍 数 。 
基于 式 (3)， 上 述 接收 信号 的 协 方差 可 表示 为 
R=E{X(i)X"() 
=APA” +W 4) 
= 及 ,+ 全 
及 为 接收 信号 协 


其 中 方差 和 矩阵， 


L 人 一 
P= 了 》 SS =diag {P,P,%, P,P,) 为 埋 号 功率 协 方差 矩 车 ， 


t=!1 


PP 为 单 信号 功率 ，R, = APA” 为 不 包含 噪声 的 信号 协 方差 矩阵 ， 


非 均 义 高 斯 噪声 功率 协 方 差 矩 阵 
厂 =E{N(DJN (=diagfaa poco ， 且 


噪声 功率 ， 信 


O702 关 … 关 02 郑 … 关 02，02 为 单 非 均匀 高 
号 与 噪声 互 不 相关 。 

式 (4) 信 号 协 方差 矩阵 模型 可 通过 以 MUSIC 方法 为 代表 的 
子 空 间 类 算法 进行 特征 空间 分 解 实现 DOA 参数 估计 。 然 而 ， 
非 均 匀 高 斯 噪声 下 ， 信 号 协 方差 矩阵 特征 空间 分 解 会 引起 信号 
子 空间 泄露 ,诸如 子 空间 类 算法 的 性 能 将 会 急剧 下 降 甚 至 失效 。 
2 ”基于 低 秩 和 矩阵 恢复 的 无 噪声 信号 协 方差 重 构 算 法 

假设 阵列 信号 模型 中 的 阵 元 数 远 大 于 信号 数 〈( M4 > 0 )， 
则 无 噪声 信号 协 方差 的 秩 rank(R,)=Q<M ， 即 无 噪声 信号 协 
方差 矩阵 丸 是 一 个 低 秩 矩阵 ， 且 协 方差 矩阵 R 非 主 对 角 线 元 
素 值 等 价 于 接收 信号 对 应 元 素 值 ， 故 可 进一步 通过 低 秩 矩阵 恢 
复 理论 对 R, 进行 主 对 角 线 元 素 重 构 , 以 此 消除 非 均匀 高 斯 噪声 
对 信号 协 方差 R 的 影响 。 


对 给 定 矩阵 Y 。 Co ， 其 在 相应 的 子 集 O 内 投影 ， 即 采样 
09 可 表示 为 
X, (i,j)eQ 
P(X)) =1 5 
[ ol 小 0 otherwise ©) 


其 中 : 子 集 Qc [m]x[n] 表示 采样 元 素 案 引 集 合 ,m,n 分 别 为 采 


样 最 大 案 引 行 数 和 列 数 ，X,,[Ps(X)], 分 别 为 矩阵 天 和 


P,() 第 六 个 索引 元 素 ，P,() 为 正 交 投影 算 子 。 

鉴于 子 集 0 是 随机 抽取 的 ,使 得 采样 后 的 信号 空域 数据 秆 
阵 满足 矩阵 低 秩 特性 ， 所 以 可 通过 低 秩 箱 阵 恢复 方法 重 构 出 只 
中 的 未 知 元 素 ， 即 重 构 模型 可 表示 为 


min rank( 脏 ) 
(0) 
st. Po(¥)=Pa(X) 
其 中 : 这 为 待 重 构 变 量 ， 盛 为 其 对 应 的 已 知 变量 。 
由 于 rank( 肚 ) 秩 函 数 的 非 凸 性 使 得 上 式 难 以 求解 。 为 求解 
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为 便于 后 续 低 秩 和 矩阵 恢复 算法 求解 ， 式 (15) 凸 优化 求解 问 
步 表示 为 


min 上 
Ro,t 


St. «(r+ 十 一 RR 区 


|Jvec(R, -R)D <¢ <E 
其 中 : 1 为 辅助 优化 变量 。 
优化 问题 (16) 中 约束 条 件 可 进 


(16) 


站 


步 改 写 为 
«(SR <1-u(rR,) 
| Jvec(R, -R)] [Jvec(R, -R)| < Ee 


(17) 


根据 Schur 补 定理 1 及 tr (Xx*x)=[vec(X)] vec(X) 可 知 ， 


式 (17) 凸 优化 求解 问题 可 转换 为 SDP 问题 ， 即 


区 


min 上 
Ro,t 


| 
党 友 


(r = tr(rRo)) (vec(R,)) 


录用 稿 
问题 式 (6)， 可 将 其 松弛 为 凸 核 范 数 | 圣 上 | ， 即 
in | 训 
mm || 
st. Pao(¥)=P,(X) 
其 中 核 范 数 | 玉 | 等 价 于 矩阵 不 所 有 奇异 值 之 和 。 
为 避免 利用 低 秩 矩 阵 恢复 算法 对 强 相关 性 数据 进行 求解 时 
出 现 解 不 稳定 的 情况 ， 本 文 引入 弹性 正则 化 项 | 
矩阵 恢复 算法 求解 稳定 性 Ho, 则 式 (7) 凸 优化 模型 进一步 表示 为 : 
mi oz|+3lzE 攻 
st. Po(¥)=Po(X) 
其 中 ， = 为 | 如 和 | 之 间 的 平衡 度 正则 化 因子 。 
基于 式 (4)(8)， 低 秩 和 矩阵 恢复 求解 模型 可 改写 为 
1 2 
min rR + le (9) 
st. Po(R,)= Po(R) 
无 噪声 协 方差 R, 是 一 个 半 正 定 和 矩阵 可 得 
rR =7tr (Ro)=tr (rR,) (10) 


其 中 : tr(.) 表示 和 矩阵 的 迹 ， 即 矩阵 对 角 线 元 素 之 和 。 


又 |R =[(RR -VC | =u(ReR,)’ RO) 


vec(R,) I (18) 
(Jvec(R, —R)) -0 
Jvec(R, 一 下) I 
上 述 SDP 问题 可 运用 CVX 工具 包 实 现 高 效 求解 上 。 
由 上 述 分 析 知 ,本 文 基于 低 秩 矩阵 恢复 理论 重 构 尺 ， 以 消 


除非 均匀 
数 估计 。 


噪声 影响 , 而 后 可 采 


传统 子 空间 类 算法 实现 DOA 参 


最 小 化 项 可 重 写 为 


rm + 


其 中 丁 代表 一 个 M(M -1)xM? 的 选择 矩阵 。 
基于 式 (11)(12)， 式 (9) 凸 优化 模型 可 


同时 ， 式 (9) 约 束 条 


Rl = (rR )+3t (RR,) 


-rR + | 


件 P,(R,)=P 
Jvec(R,—R)=0 


| “(en, IRR, 


式 (13) 凸 优化 模型 


st. Jvec (R, 一 玉 ) =0 


型 是 一 个 线性 规划 问题 ， 约 


进一步 改写 为 


(11) 


,(R) 可 等 价 表示 为 


(12) 


束 条 件 


程 中 要 求 较为 苛刻 故 可 [将 


其 等 价 松 


Jvec(R, 一 玉 ) =0 在 实际 


弛 化 为 


其 中 : & 表示 为 与 接收 信 


Jvec(R, — RR) <¢ 


则 式 (14) 凸 优化 模型 可 等 价 表示 为 : 


min oh, rasR, | 


St. |mec(R。 -RD <é 


(14) 


号 协 方差 R 相关 的 一 误差 常数 。 


(15) 


然而 ， 需 要 注意 的 是 ， 基 于 R 的 传统 DOA 估计 算法 
没有 考虑 到 协 方 差 矩阵 R 元 素 之 间 的 的 相关 性 ， 在 实现 DOA 
参数 估计 的 同时 ， 会 导致 算法 计算 量 较 大 、 低 信 品 比 条 件 下 估 
计 稳定 性 差 和 分 辩 率 低 等 问题 。 


3 ”二 阶 统计 量 域 下 的 加 权 上 L1 稀疏 重 构 算法 


为 进一步 提高 低 信 噪 比 条 件 下 信和 号 估计 精度 和 分 辨 为 ， 降 
低 稀 疏 重 构 算法 计算 复杂 度 ， 本 文 基于 低 秩 矩 阵 恢复 理论 ， 提 
出 一 种 二 阶 统计 量 域 下 加 权 7 范 数 稀疏 重 构 算 法 。 

无 噪声 协 方差 R, 中 的 元 素 ， 即 第 ; 行 和 第 j 列 的 阵列 输出 
互相 关系 数 可 表示 为 


广 =(R 0); =(APA”), -Dae 站 dsingl /4 (19) 
对 信号 协 方差 R, 按 行 展 开 ， 可 表示 为 
ni Tm 
R -| : : 
Tal Tuam 
Ee j2x(M -1)d sin, /4 
和 Be Ye | (20) 


g=1 g=1 


Q 
—j2z(M -1)dsinQ, /4 
2 Pe 


ql 4=1 


其 中 : 当 i- j=u-v 时 ， 互 相关 系数 = ，i,j,uve[l,M]; 
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即 它们 含有 相同 系数 特征 ， 可 通过 二 阶 统计 量 域 下 求 和 平均 将 
多 矢量 问题 转换 为 单 矢量 问题 。 

定义 (2M -TJxl 的 向 量 Y ， 通 过 对 信号 协 方差 R, 求 和 平 
均 ， 其 第 而 个 元 素 可 表示 为 


YR (iM +n 一 万 )， m=1,…,M 
了 (元 )= ee Q21) 
1 2M -页 
PE > R(xrmy),m=M +L…2M -1 
一 171 j=1 
其 中 ，mr=M+1- 元 六 =2M +1- 元 一 而 。 
了 的 向 量 形式 可 表示 为 
r=B(0)P (22) 
其 中 : 已 =|B, 忆 已 | ， 
B(0)=[5(0),…,b(0,),…b(0,) |eC™ 7 
为 虚拟 阵列 流 型 矩阵 ， 其 第 gq 列 可 表示 为 
bp(9 = hn A ns 
0)-[ 2 


C127(M 2)dsin6, /4 20600 
， ， 


区 


单 矢 量 问题 可 转变 为 稀 疏 信号 重 构 问 题 ， 将 空间 波 达 


方向 角 9 分 割 为 V (N > M) 个 网 格 , 即 到 = 仿 ,@ ,6 风 


稀 琉 后 的 接收 信号 矢量 可 表示 为 
六 = 万 (到 ) 石 (24) 


其 中 ，B()=|5(6),b(&)…,5(B,) |e CW* 为 过 完备 基 抵 


阵 ，=[B, 避 ,…,B | 为 稀疏 的 信号 功率 向 量 。 


上 述 问 题 可 通过 7 范 数 约束 最 优化 算法 进行 求解 ， 即 


min|B| si |r_B(P)B), <7 05) 


其 中 : 为 误差 参数 因子 ， ,为 2 范 数 和 1 范 数 。 
业已 证 明 : 优化 问题 (25) 中 7 范 数 约束 对 目标 信号 系数 呈 

现 不 均匀 惩罚 ， 其 对 大 系数 惩罚 程度 大 于 小 系数 ， 导 致 其 事实 

上 是 一 种 有 偏 的 估计 031。 

针对 此 问题 ， 本 文采 用 加 权 7 范 数 重 构 稀 疏 信号 功率 谱 ， 

而 后 通过 谱 峰 搜 索 寻 找 最 大 谱 峰 对 应 信 源 位 置 从 而 实现 非 均匀 

噪声 下 DOA 估计 ， 即 


F 和 


min|w oP] srr -B(P)F| 


2 <7 (26) 
式 (26) 可 进一步 改写 为 
min Sw, 


式 (27) 凸 优化 问题 可 通过 二 阶 锥 规划 (Second-order cone 
programming, SOCP ) 04 重 构 实 现 DOA 参数 估计 。 
基于 上 述 讨 论 ， 本 文 提出 的 LR-WLOSRSS 算法 可 表述 如 


Blst. |r-B(Y)B| <7 (27) 


下 : 
初始 化 
1l. Y=avr(Ro), Maxiter, €£», XX», 71; 


2 Buemin FAlst pr -BCP)B), <7 

迭代 | 

3. 权 值 w=J/(B (i)+X): 

4 求解 如 下 优化 间 是 


B,,= a Yi | 下 sz |Y-B()Bl < 以 更 新 系数 ; 
i=l 


5. 重复 步骤 3、4 直至 满足 如 下 条 件 | ,~ 户 , ,| < 或 者 


I 


m2> Maxiter 。 

其 中 算 子 avr(.) 代表 式 21) 矩 阵 求 和 平均 ，Maxiter 表示 最 大 
友 代 次 数 ， s 为 迭代 终止 参数 ，%% 为 权 值 w 平衡 约束 因子 ， 防 
止 算法 迭代 时 权 值 分 母 趋 于 无 穷 大 ; 7 为 凸 优化 问题 误差 参数 ， 


m 为 迭代 次 数 ，B, ,为 谱 峰 估计 初始 值 ; ,为 第 ma 次 迭代 的 


谱 峰 估计 值 。 
4 ”实验 仿真 及 分 析 


本 节 通 过 实验 仿真 分 析 进 一 步 验证 所 提 算 法 的 有 效 性 。 实 
仿 仿 真 条件 如 下 : 均匀 线性 阵列 空间 阵 元 个 数 M =8， 信 号 快 
拍 数 L=500 ， 非 相干 信号 功率 为 o? ， 误 差 常数 &=5 ，z=60， 
a=0.01， 和 =0.1P, ， 网 格 数量 N=181 。 信 噪 比 0 定义 为 


M 2 
SNR= 一 名 ， 其 中 为 2 单 噪声 功率 。 均 方 根 误差 定义 为 
ml Om 
(5 -CT 08) 
RMSE = | 一 一 = 
3 gk ) 


其 中 : 为 蒙特 卡 洛 实验 次 数 。 

非 均 匀 高 斯 噪声 功率 协 方差 定义 为 

W=diag{2.0,10,2.5,5.0,0.5,1.5,3.0,5.0} (29) 

其 中 : diag {®} 表示 对 角 化 算 子 。 

实验 1 考虑 SNR=0dB 和 -5dB 两 种 情况 。 图 1 (a) 为 入 
射 角度 分 别 为 -3*、10°* 和 16°%, SNR=0dB 的 空域 谱 。 由 图 1 (a) 
知 ，MUSIC、WL1 和 1-SVD 算法 在 非 均匀 噪声 条 件 下 不 能 
确 分 辨 位 于 10* 和 “16° 的 两 个 目标 ， 而 MC-MUSIC 和 LR- 
WLOSRSS 算法 可 正确 分 辨 10" 和 16? 的 目标 。 图 1 (b) 为 入 射 
角度 分 别 为 -3"、10* 和 16", 且 SNR=-5dB 的 空域 谱 。 由 图 1 (b) 
知 ， 低 信 噪 比 情况 下 ， 本 文 所 提 LR-WLOSRSS 算法 可 有 效 分 
辨 三 个 角度 ， 而 MC-MUSIC 等 其 他 算法 均 不 能 分 辨 10" 和 16。 
两 个 目标 。 图 1 表明 , 与 MC-MUSIC 等 其 他 传统 类 算法 相 比 ， 
所 提 算 法 在 非 均匀 噪声 和 低 信 噪 比 条 件 下 仍 具 有 较 好 的 DOA 
估计 性 能 。 另外 , 从 图 1 还 可 看 出 , 本文 所 提 LR-WLOSRSS 算 
法 具有 较 窄 的 主因 及 较 低 的 旁 六 ， 表 明 所 提 算 法 在 非 均 匀 高 斯 
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噪声 和 低 信 品 比 条 件 下 具有 较 高 的 角度 估计 精度 和 分 辨 力 。 


MC-MUSIC 


— MC-MUSIC 


MUSIC 
LR-WLOSRSS, 
1 


LISVD 


归 一 化 输出 功率 /dB 
归 一 化 输出 功率 /dB 


(a) SNR=0dB 


图 1 信 噪 比 SNR=0dB 和 -5dB 条 件 下 的 非 相干 信号 空域 谱 对 比 图 


(b) SNR=-5dB 


实验 2 考虑 入 射 角度 分 别 为 -3°*、10°* 和 13° 的 非 相 干 信号 ， 
SNR=5dB。 图 2 为 五 种 算法 空域 谱 。 由 图 2 可 知 ，MUSIC 和 
WL1 算法 不 能 有 效 分 辨 10" 和 13° 两 个 邻近 角度 ， MC-MUSIC 
虽 可 重 构 出 无 噪声 协 方差 ， 但 受 限 于 算法 角度 分 辩 率 影响 ， 不 


房 云 飞 ， 


归 一 化 输出 功率 /dB 


-60 -40 20 0 20 40 60 -60 40 0 
DOA/(®) DOA/(°) 
(a) (b) 


图 3 信 源 数 变化 对 MUSIC、LR-WLOSRSS 算法 的 影响 


实验 4 考虑 两 个 入 射 角度 为 -3* 和 5° 的 非 相 干 信号 , 快 拍 
数 L=500， 信 噪 比 SNR=[-8:2:12]， 进 行 200 次 蒙特 卡 洛 独立 重 
复 实验 。 由 图 4 和 表 2 可 知 ， 传 统 MUSIC 算法 在 低 信 噪 比 条 
件 下 的 DOA 估计 RMSE 相对 较 高 。MC-MUSIC 算法 采用 矩阵 
补 全 理论 重 构 出 无 噪声 协 方差 ， 以 消除 非 均匀 噪声 的 影响 ， 故 


能 对 邻近 角度 实现 有 效 估计 。 本 文 所 提 算 法 及 理想 高 斯 白 噪 
下 的 1-SVD (Ideal-1-SVD 〉 可 对 邻近 角度 实现 有 效 分 辨 。 需 
要 注意 的 是 ，Jdeal-1-SVD 虽 可 正确 估计 出 两 邻近 角度 , 但 相对 
于 LR-WLOSRSS 而 言 ，1deal-1-SVD 在 两 邻近 角度 之 间 的 谱 峰 
估计 精度 较 差 ， 而 LR-WLOSRSS 则 估计 精度 较 高 ， 且 具有 更 
低 旁 汶 , 表 明 LR-WLOSRSS 对 邻近 角度 具有 较 高 角度 分 辨 力 。 


MC-MUSIC 


USIC 
LR-WLOSRSS 
WL1 


ldeaH-L1-SVD 


归 一 化 输出 功率 /dB 


60 -40 -20 20 40 60 


0 
DOAI(9) 


图 2 邻近 非 相 干 信号 空域 谱 对 比 图 


实验 3 考虑 三 个 入 射 角度 分 别 为 -3*>、10°* 和 16° 的 非 相 干 
信号 ， 信 噪 比 SNR=5dB， 信 源 数 K=1，2，3 时 , MUSIC、LR- 
WLOSRSS 算法 的 DOA 估计 有 效 性 。 图 3 (a) 为 信 源 数 变 化 
对 MUSIC 算法 性 能 影响 。 从 3 (a) 可 以 看 出 ， 当 信 源 数 K 分 
别 为 1，2 时 ， 传 统 MUSIC 算法 不 能 有 效 分 辨 三 个 角度 。 只 有 
当 信 源 数 等 于 待 估计 角度 数量 ， 即 K=3 时 ，MUSIC 才能 有 效 
分 辨 三 个 角度 。 图 3 (b) 为 信 源 数 变 化 对 LR-WLOSRSS 算法 
性 能 影响 。 从 3 (b) 可 以 看 出 , 当 信 源 数 K 分 别 为 1, 2, 3 时 ， 
本 文 所 提 LR-WLOSRSS 算法 均 能 有 效 分 辨 三 个 角度 ， 信 源 数 
的 变化 对 本 文 所 提 LR-WLOSRSS 算法 几乎 没有 影响 。 从 图 3 
可 以 看 出 ， 传 统 子 空间 类 DOA 估计 算法 由 于 受制 于 信号 或 信 
号 子 空间 数目 ， 即 信 源 数目 的 先 验 信息 ， 导 致 子 空间 类 算法 的 
DOA 估计 性 能 对 信 源 数 的 变化 较为 敏感 ;本文 所 提 LR- 
WLOSRSS 算法 利用 目标 空域 稀疏 性 ， 采 用 稀 玻 重 构 对 待 估计 
信号 进行 稀 玻 逼近， 而 无 须 已 知 信 源 先 验 信 息 。 


oe 


而 MC-MUSIC 算法 在 低 信 噪 比 条 件 下 DOA 估计 RMSE 低 于 
传统 MUSIC 算法 。 在 低 信 噪 比 条 件 下 ，WL1 算法 和 1-SVD 算 
法 有 一 个 相对 较 低 的 RMSE。 此 外 , 由 图 4 和 表 2 还 可 以 看 出 ， 
在 给 定 仿真 条 件 下 ， 本 文 所 提 LR-WLOSRSS 算法 估计 性 能 
RMSE 意义 下 明显 优 于 WL1 算法 和 1-SVD 算法 ， 在 低 信 噪 比 
条 件 下 ， 此 优势 尤为 明显 ， 表 明 所 提 LR-WLOSRSS 算法 具 
较 好 的 DOA 估计 性 能 。 


—+— MUsIC 

一 举 一 MC-MUSIC 
二 

一 6 一 LR-WLOSRSS 
一 qd 一 LI-SVD 


in 


估计 均 方 根 误差 /(*) 


4 1 
8 -6 -4 2 


0 2 4 6 8 10 12 
信 噪 比 (SNR) /dB 


图 4 DOA 估计 RMSE 随 信 噪 比 SNR 的 变化 对 比 图 


表 1 DOA 估计 RMSE 随 信 噪 比 SNR 变化 


》 SNR/dB 
算法 -8 =6 -4 0 4 8 12 
MUSIC 2.4098 2.3965 2.3360 1.0588 0.2125 0.1010 0.0650 


MC-MUSIC 1.8698 1.1065 0.8060 0.3009 
WL1 1.4698 0.7495 0.4760 0.1509 
LR-WLOSRSS 0.9098 0.4695 0.2060 0.1085 
LI-SVYD _1.1698 0.7265 0.4260 0.1709 


0.1625 0.1007 0.0641 
0.0726 0.0625 0.0601 
0.0725 0.0624 0.0600 
0.0725 0.0625 0.0600 

实验 5 考虑 两 个 入 射 角度 为 -3° 和 5? 的 非 相干 信号 ， 信 噪 
比 SNR=0dB， 快 拍 数 L=[100:1200]， 进 行 200 次 蒙特 卡 洛 独立 
重复 实验 。 由 图 5 和 表 2 可 知 ， 随 着 快 拍 数 的 增加 ， 本 文 所 提 
LR-WLOSRSS 算法 及 MUSIC、MC-MUSIC、WL1 和 1-SVD 算 
法 的 RMSE 逐渐 降低 。 然 而 ,需要 注意 的 是 ,所 提 LR-WLOSRSS 
算法 RMSE 均 低 于 MUSIC、MC-MUSIC、WL1 和 1-SVD 算法 ， 
表明 所 提 LR-WLOSRSS 算法 在 非 均匀 高 斯 噪声 下 具有 更 优 的 
DOA 参数 估计 性 能 。 
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3 - 
-一 一 MUSIC 
一 举 一 MC-MUSIC 
2.5 一 全 一 WLI 
一 6 一 LR-WLOSRSS 
一 生 一 LI-SVD 
9 200 300 400 500 600 700 800 900 1000 1100 1200 
快 拍 数 
图 5 DOA 估计 RMSE 随 快 拍 数 的 变化 对 比 图 
表 2 DOA 估计 RMSE 随 快 拍 数 变化 
算法 快 拍 数 
、 100 300 500 700 900 1100 1200 
MUSIC 2.4398 2.1098 1.0588 0.7088 0.6059 0.5359 0.5259 
MC-MUSIC 1.9698 0.8106 0.3509 0.2125 0.1401 0.1301 0.1141 
WL1 1.7698 0.6926 0.2509 0.1429 0.1193 0.1046 0.0991 
LR-WLOSRSS 0.9698 0.4560 0.1645 0.1355 0.1150 0.0971 0.0871 
LI-SVD 1.4698 0.5560 0.2409 0.1425 0.1175 0.1011 0.0911 
土 十 
5 ”结束 语 


为 提高 非 均匀 噪声 下 DOA 估计 算法 的 估计 精度 和 分 辩 率 ， 
本 文 基于 低 秩 矩阵 恢复 理论 ， 提 出 了 一 种 三 阶 统 计量 域 下 的 加 
权 7 稀疏 重 构 DOA 估计 算法 。 该 算法 基于 低 秩 矩阵 恢复 理论 ， 
引入 弹性 正则 化 因子 将 接收 信号 协 方差 矩阵 重 构 问 题 转换 为 可 
获得 高 效 求解 的 SDP 问题 , 而 后 在 二 阶 统 计量 域 下 通过 算 阵 求 
和 平均 将 稀疏 多 矢量 重 构 问 题 转换 为 单 矢 量 重 构 问 题 ， 最 后 采 
用 加 权 刀 范 数 对 稀 疏 单 矢量 进行 重 构 以 实现 DOA 参数 估计 。 
仿真 结果 表明 , 与 传统 的 MUSIC、MC-MUSIC、WL1 和 1-SVD 
算法 相 比 ， 所 提 算 法 在 非 均 匀 高 斯 噪声 和 低 信 品 比 条 件 下 ， 具 
有 较 好 的 DOA 参数 估计 性 能 。 


参考 文献 : 


[1] Van Trees, Harry L. Optimum array processing: part IV of detection, 
estimation and modulation theory [M]. Hoboken: Wiley, 2002. 
[2] Schmidt R. Multiple emitter location and signal parameter estimation [J]. 


IEEE Trans on Antennas and Propagation, 1986, 34 (3): 276-280. 


. hinaXiv 合 作 期 十 
房 云 飞 ， 等 : 基于 低 秩 和 矩阵 恢复 的 DO 厅 ; 


[3] Liao B, Huang L, Guo C，et al. New approaches to direction-of-arrival 
estimation with sensor arrays in unknown nonuniform noise [J]. IEEE 
Sensors Journal, 2016, 16 (24): 8982-8989. 

[4] Malioutov D, Cetin M, Willsky A S$. A sparse signal reconstruction 
perspective for source localization with sensor arrays [J]. IEEE Trans on 
Signal Processing, 2005, 53 (8): 3010-3022. 

[5] Pesavento M, Gershman A B. Maximum-likelihood direction-of-arrival 
estimation in the presence of unknown nonuniform noise [J]. IEEE Trans on 
Signal Processing, 2001, 49 (7): 1310-1324. 

[6] Liao B, Guo C, Huang L, et al. Matrix completion based direction-of-arrival 
estimation in nonuniform noise [C1]// Proc of IEEE International Conference 
on Digital Signal Processing. 2016: 66-69. 

[7] Candes EJ, Recht B. Exact matrix completion via convex optimization [J]. 
Foundations of Computational Mathematics, 2009, 9 (6): 717. 

[8] Candes E J, Plan Y. Matrix completion with noise [J]. Proceedings of the 
IEEE, 2010, 98 (6): 925-936. 

[9] Candes E J, Wakin M B, Boyd S P. Enhancing sparsity by reweighted L1 
minimization [J]. Journal of Fourier Analysis and Applications, 2008, 14 (5): 
877-905. 

[10] Cai J F, Candes E J, Shen Z. A singular value thresholding algorithm for 
matrix completion [J]. SIAM Journal on Optimization, 2010, 20 (4): 1956- 
1982. 

[11] Horm R A, Johnson C R. Matix analysis [M]. Cambridge: Cambridge 
university press, 2012. 

[12] Grant M, Boyd S$, Ye Y. CVX: Matlab software for disciplined convex 
programming, version 2. 1 [EB/OL]. [2017-03-01]. http://cvxr. com/cvx. 

[13] Tian Y, Sun X, Zhao S. DOA and power estimation using a Sparse 
representation of second-order statistics vector and LO-norm approximation 
[J]. Signal Processing, 2014, 105: 98-108. 

[14] Lobo M S, Vandenberghe L, Boyd S, et al. Application of second-order cone 
programming [J]. Linear Algebra and its Applications, 1998, 284 (1-3): 193- 
228. 

[15] Liao B, Chan S C, Huang L, et al. Iterative methods for subspace and DOA 
estimation in nonuniform Noise [J]. IEEE Trans on Signal Processing, 2016, 


64 (12): 3008-3020. 


